Compressing Random Microstructures via Stochastic Wang Tilings 
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This paper presents a stochastic Wang tiling based technique to compress or reconstruct disordered 
microstructures on the basis of given spatial statistics. Unlike the existing approaches based on a 
single unit cell, it utilizes a finite set of tiles assembled by a stochastic tiling algorithm, thereby 
allowing to accurately reproduce long-range orientation orders in a computationally efficient manner. 
Although the basic features of the method are demonstrated for a two-dimensional particulate 
suspension, the present framework is fully extensible to generic multi-dimensional media. 
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In 1961, Hao Wang introduced a tiling concept based 
on square dominoes with colored edges permitting their 
mutual assembly in a geometrically compatible (hard) 
manner pQ. Since then, his tiles have been the subject of 
studies in discrete mathematics [2H4] and found an ex- 
tensive use in computer graphics [5 , game industry [6], 
theory of quasicrystals [7 or biology [8]. From the per- 
spective of this paper, the appealing feature of Wang 
tilings is that they can compress and reproduce naturally 
looking planar patterns or three-dimensional surfaces by 
employing only a small number of distinct tiles [3j [9j [10] . 
Motivated by this observation, we further explore the po- 
tential of Wang tiles to represent long-range spatial cor- 
relations in disordered microstructures; a problem com- 
mon to materials science [IT], geostatistics [12] or image 
analysis [15] . 

In this regard, two closely related applications can 
be distinguished, namely the micro structure reconstruc- 
tion [T4HT6] based on given spatial statistics and mi- 
cro structure compression jT7HT9] aiming at efficient rep- 
resentation of materials structure in multi-scale compu- 
tations [20]. Our focus is on the latter, since these pro- 
cedures usually have the microstructure reconstruction 
techniques at heart, hereby covering the common fea- 
tures of both. 

To the best of our knowledge, compression algorithms 
reported to date use a single cell (PUC) that is peri- 
odically extended to tile the plane in a deterministic 
manner [18]. Such structures then inevitably manifest 
strong long-range correlations with a period of the PUC 
dimensions. We shall demonstrate that these artifacts 
can be effectively controlled when utilizing small Wang 
tile sets [5J [10] , carefully designed to capture morphologi- 
cal trademarks of compressed media, combined with fast 
stochastic Cohen-Shade-Hiller-Deussen (CSHD) tiling al- 
gorithm [5 for real-time texture generation. A poten- 
tial of this approach will be demonstrated for equilib- 
rium two-dimensional particulate suspensions consisting 



of equi-sized disks of radius p uniformly distributed in a 
homogeneous matrix, cf. [2T] , 

To this goal, consider a two-dimensional microstruc- 
tured domain V discretized by a regular square lattice. 
Each lattice cell contains specific morphological patterns 
that are compatible on contiguous boundaries, Fig. [TJb). 
If there are no missing cells in V, the discretization is 
called a valid tiling, and a single cell is referred to as the 
Wang tile [1 , Fig. [TJa). The tiles have different codes 
on their edges, lower-case Greek symbols in Fig. [TJa), 
and are not allowed to rotate when tiling a plane. The 
number of distinct tiles within V is limited, though ar- 
ranged in such a fashion that none of them or any of their 
sub-sequence periodically repeats. The gathered distinct 
tiles are referred to as the tile set, Fig. [TJa). Sets that 
enable uncountably many, always aperiodic, tilings are 
called aperiodic [5JH]. In real world applications, the as- 
sumption of strict aperiodicity of the tile sets is relaxed 
to aperiodicity of tilings, ensured e.g. by the CSHD al- 
gorithm [5] introduced next. 
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FIG. 1. (a) A Wang tile set W8/2-2 with edge length £ and 
codes {a, /3,7, S}. (b) Example of an aperiodic valid tiling. 

Intuitively, the ability of a tile set to control long-range 
order effects arises from tile and edge code diversities, 
Fig. []Ja). Both factors are related, so that while the 
number of edge codes n\ in z-th spatial direction can be 
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chosen arbitrarily, the number of tiles n 1 must satisfy 

n l = n NW Vn^, (1) 

where n cs = (n^n^) 2 is the number of tiles in the com- 
plete set and n NW = 2, . . . , \fn^ stands for the optional 
number of tiles with identical arrangement of north- 
western (NW) edge codes. The complete set of n cs tiles 
is obtained by permuting the chosen codes c$. In valid 
tilings, the south-eastern edge codes must match those 
assigned to NW edges, Fig. [ljb) . Thus, the tiles in the 
complete set are sorted according to existing NW edge 
code combinations and the desired number of tiles in a 
user defined set is formed by selecting n NW of unique tiles 
from each group. Such a set is referred to as Wn* /n^-n^ 
in what follows. Notice that the Wl/1-1 set corresponds 
to the PUC setting. 

In the stochastic tiling algorithm, the index of a new 
tile to be placed is selected randomly with the uniform 
probability from an appropriate NW group compatible 
with the eastern code of the tile previously placed and the 
southern code of the tile above the one to be placed (edges 
7 and a in bold adjacent to shaded cell in Fig. [TJb)). 
Aperiodicity of the resulting tiling is guaranteed provided 
that the random generator never returns a periodic se- 
quence of numbers and that each NW group contains at 
least two distinct tiles [5]. 

Analogously to the existing works on reconstruction 
and compression of random media, the tile morphologies 
are designed by an optimization procedure expressed in 
terms of suitable statistical descriptors. As our focus is 
to control long-range artifacts, we limit the exposition to 
the two-point probability function S 2 (x) [TT]. For statis- 
tically uniform ergodic media, it provides the probability 
that two arbitrary points from X>, separated by x, are 
simultaneously found in the particle phase. The function 
satisfies 62(C)) = </>, where (j) is the volume fraction of 
particles, and S2 (x) « (j) 2 for ||x|| > A indicates the ab- 
sence of long-range orders at the characteristic length A, 
Fig. lb). 

In the current setting, the Wang tiling compression 
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FIG. 2. (a) Reference two-phase medium of size 174p x 174p 
formed by equilibrium distribution of 1,300 equi-sized disks 
of volume fraction 26.8% and (b) its two-point probability 
function S2] p is the disk radius. 



consists of a set of n 1 tiles of the edge length £, in which 
we distribute n d disks. The configuration of particles 
is determined by the parameter vector £2,^=1? 
where td G {1, ... ,71*} specifies the parent tile index of 
the d-th disk and £j jC j G [0, £] the local position of the disk 
at j-th direction. To determine the two-point probability 
function S2 for a given configuration, we assemble a tiling 
that covers the domain of the same size as the represen- 
tative sample P, Fig. |2|a). Notice that such tiling cor- 
responds to a realization of a statistically homogeneous 
material, since the tiles are selected from NW edge groups 
with the uniform probability. The proximity of the tile- 
based morphology to the original sample is quantified by 
an objective function 

£=p^(s 2 (x)-5 2 (x)) 2 dx (2) 

which can be efficiently evaluated using the Fast Fourier 
Transform techniques, e.g. [11 . The minimization of ([2| 
is carried out by the Simulated Re-Annealing method 
with computational cost similar to existing PUC design 
strategies [22 . The algorithm ensures that the tiles in 
the set satisfy the corner constraint [5] , requiring that the 
tile corners are not occupied by a disk, and determines 
the number of disks n d and the cell size i such that the 
local volume fractions associated with edges (grey disks 
in Fig. |5| and tile interiors (white disks in Fig. [5| are as 
close to the target value <p as possible. 

An example of optimal approximations of the target 
microstructure from Fig. [2] in terms of a PUC and the 
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FIG. 3. Optimized microstructures and two point probability 
functions S2 for PUC (a,b) and set W18/3-3 (c,d). The arrow 
in (c) denotes periodic region due to local character of tile 
placement in CSHD tiling algorithm. 
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FIG. 4. (Color online) Comparison of & in X1S2 plane for (a) tile set W8/2-2 with respect to number of disks n d involved, 
(b) different sets, (c) Comparison of different sets in terms of objective function (pi) and with respect to number of disks n d 
involved. The curves in (b) are plotted for particular values of n d highlighted in (c)7 



Wang tile set W18/3-3 is shown in Fig. [3J The repre- 
sentations are based on n d = 49 particles and tile sizes 
i = 24.5 p and t = 7.5 p, respectively. Evidently, both 
heterogeneity patterns carry long-range order effects with 
the period of i, manifested as the local peaks 

S 2 = maxS 2 (k£) (3) 

k\{0} 

in the two-point probability functions, Fig. [3jb,d). No- 
tice that S2 is always equal to <\> for the PUC approach, 
whereas the Wang tiles are capable of adjusting these 
artifacts by the proper morphology design. This is also 
reflected in visual regularity of the generated suspensions, 
compare Fig. ga) with Fig. [3la,c). Also observe the lo- 
cally periodic region in Fig. Mc), arising from the local 
character of CSHD algorithm and from the lowest num- 
ber of tiles in groups of admissible NW edge code com- 
binations, n NW = 2. Such phenomenon is thus less likely 
when increasing parameter n NW , however, at the expense 
of increasing set sizes, especially for higher edge code di- 
versities, recall Eq. ([!]). 

(a) (b) 

FIG. 5. Building blocks of microstructure compression based 
on (a) PUC and (b) tile set W18/3-3 with 49 disks assigned 
to tile edges (grey) and interiors (white). 

The principal features of Wang tile-based compressions 
are further illustrated in Fig. [5] Instead of relying on a 
single cell containing the complete morphological infor- 
mation, Fig. ga), the tiling-based approach utilizes sub- 
stantially simpler building blocks, Fig. [5|b), assembled 
to comply with edge constraints (grey disks in Fig. [5|. 
This, however, restricts the space of admissible disk con- 
figurations in Wang tiles compared to the single PUC 
design. 



Now, we are in the position to quantify to which ex- 
tent is the quality of reconstructed suspensions deter- 
mined by the tile set diversity and the morphology de- 
sign itself. This aspect is examined first in Fig. [4] by 
means of sections of the two-point probabilities 62(^1, 0) 
and the objective function E, revealing that two effects 
govern the amplitude and period of the local extremes 
S 2 - First, for a fixed tile set, increasing the number of 
disks increases the tile edge dimension t (and thus the pe- 
riod) and slightly decreases the amplitude, Fig.ga). On 
the other hand, increasing the number of tiles decreases 
the period as well as the magnitude of local extremes, 
Fig. gb). Also notice that the quality of the W18/3-3 
set in terms of the objective function ^ is systemati- 
cally inferior to W8/2-2, Fig. gc). This is caused by 
an inaccurate representation of disk volume fraction for 
the former set, which pollutes the shape of 62 statistics, 
Fig. gb) [23 . It further follows from Fig. gc) that in- 
creasing the tile set diversity is much more efficient; for 
sets containing more than 32 tiles, the error is almost in- 
dependent of the number of disks. This saturation value 
reflects the inaccurate representation of short-range val- 
ues of S2, caused by the particular form of the objective 
function ([2]). If needed, the local details can be incor- 
porated in terms of higher-order statistics or specifically 
tailored descriptors [T4l [T5l [241 125] . 

It is now clear that the local extremes can be attributed 
to a limited number of tiles used in a repetitive, although 
random fashion. Actually, two components repeat when 
tiling the plane: tile edges and interiors. To study the 
local artifacts analytically, we consider user defined sets 
with tiles selected so that their edges incorporate each 
code Ci at least once. Assuming that tiles and edges 
repeat independently, the maximum local extremes can 
be estimated as 

% - J [4> + (»< - 1) <?} + max { g [4> + « - 1)4?] } 
where ft = (£ — 4p) 2 /£ 2 and ft = 1 - ft denote the 
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volume fractions of tile interiors (occupied by white disks 
in Fig. [5| and edges (occupied by grey disks in Fig. [5|, 
respectively [26]. 

In Fig. |6j we compare the actual values of S2 with the- 
oretical predictions Q for several values of (f) e . Apart 
from the limit cases, (jf G {0, 1}, S% was also explored 
for (jf = 0.2 (average value from all considered tile sets). 
We observe that an almost exact match is obtained for 
the lower bound with <j) e = 0, red solid curve in Fig.|6ja), 
demonstrating that the long-range artifacts are carried 
mainly by the tile interiors. This is rather surprising, 
since all considered tile sets satisfy n 1 ^> n\ = n§ , so that 
edges repeat more often than the tile interiors. Moreover, 
the magnitude of spatial artifacts converges rapidly to the 
limit value (j) 2 . Altogether, this leads us to the conclu- 
sion that artifacts due to discrete nature of Wang tilings 
can be almost eliminated by a proper morphology opti- 
mization. Also note that the accuracy of the estimate Q 
appears to be reasonable, both for the average value of 
(jf = 0.2, blue double dotted curve in Fig.|6ja), as well as 
for values corresponding to individual tile sets, Fig.[6jb). 
It may thus serve as a basis for the a-priori selection of 
the tile set parameters n\ and n 1 . 




FIG. 6. (Color online) (a) Dependence of local extremes S2 on 
the particular tile set. (b) Correlation of local peaks obtained 
from two-point probabilities of optimized microstructures and 
their predictions given by Eq. Q; r in (b) denotes the Pearson 
correlation coefficient. 

Summary. A new compression/reconstruction tech- 
nique based on Wang tilings has been proposed and ap- 
plied to two-dimensional microstructures of disordered 
particulate media. The technique is extensible to generic 
three-dimensional microstructures, adopting the frame- 
works of Wang cubes [27j [28] and image synthesis [5] ; it 
substantially generalizes the periodic unit cell concept by 
making use of multiple tiles instead of a single cell, pre- 
serves long range spatial features, and is computationally 
efficient. A formula for estimates of long-range order spa- 
tial artifacts has also been proposed and verified for the 
studied material system. 
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